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Fundamental phase-shift detection properties of optical 
multimode interferometers are analyzed. Limits on perfectly 
distinguishable phase shifts are derived for general quantum 
states of a given average energy. In contrast to earlier work, 
the limits are found to be independent of the number of in- 
terfering modes. However, the reported bounds are consistent 
with the Heisenberg limit. A short discussion on the concept 
of well-defined relative phase is also included. 

PACS numbers: 42.50.-p, 03.65.Ta, 07.60. Ly, 42.50.Dv 



I. INTRODUCTION 

It is today well-known that the use of quantum- 
mechanical states can improve the precision of interfer- 
ometric measurements. According to the so-called stan- 
dard quantum limit Q , the precision of optical measure- 
ments employing classical states cannot increase faster 
than oc 1 /Ve, where E is the energy used in the measure- 
ment. However, the use of nonclassical states allows us to 
reach the Heisenberg limit oc 1/E which for high en- 
ergies would give a remarkable improvement in accuracy. 
Several setups have theoretically been shown to work 
at the Heisenberg limit ^^l^. The standard quantum 
limit has also been circumvented experimentally [0-0. 
However, the fragile nature of the quantum states has so 
far prevented these measurements from being carried out 
with higher energies. Therefore, high intensity classical 
interferometry reach a much better overall accuracy. 

Recently, a bound closely related to the Heisenberg 
limit was given by Margolus and Levitin . The bound 
gives the time necessary for a state of a closed system to 
become orthogonal for a given average energy. This limits 
the rate of operations in quantum information processing 
and the resolution in interferometry. In an earlier paper 
, we derived the states that minimize the time needed 
to freely evolve into another state, whose overlap with 
the original one was given. This would correspond to 
minimizing the necessary phase shift for a single-mode 
state. 

The present paper is devoted to the phase-shift detec- 
tion properties of multimode interferometers. In many 
interferometric measurements, a single induced phase 
shift is tracked by monitoring the interference fringes in 
the output of a two-mode interferometer. However, it is 
often possible to induce several phase shifts of the same, 



or smaller, magnitude (possibly with different signs) in 
the different arms of a multimode interferometer. Here, 
we investigate whether this fact can be used to improve 
the accuracy of such measurements. We will assume that 
the interfering fields have the same optical frequency, so 
that the total energy is proportional to the number of 
photons used. 

In an earlier investigation it was concluded that 
the accuracy of the multimode interferometer would im- 
prove indefinitely with the number of modes. However, 
we find that there is no fundamental advantage in using 
more than two modes and that, in our eyes, the con- 
clusions in Ref. |13| stem from an unfortunate choice of 
figure of merit. Rather, the accuracy is found to be lim- 
ited only by the energy used in the measurement and 
scale according to the Heisenberg limit. 



II. PROBLEM FORMULATION 

Accuracy can be defined in many different ways. For 
example, in a recent paper on multimode interferometry 
[Q , the width of the major peak of the phase distribution 
was taken as a measure of the precision. Here, we will 
define it as the smallest phase shift required to give a 
perfectly distinguishable outcome, i.e., the phase-shifted 
state is required to be orthogonal to the original state. 

More precisely, we look for the smallest phase shift 
that can be detected with certainty using an M-mode 
interferometer and a state whose average energy is given. 
The smallest phase shift (j) should here be interpreted in 
the following sense (cf. Figs. ^ and H). We consider the 
case where all arms of the interferometer have induced 
phase shifts that satisfy (pm = \n4'j where — 1 < Am < 
1 and (/) > is the parameter to be minimized under 
the condition that it results in a state that is perfectly 
distinguishable from the original. Since the A^-photon 
state (|0,Ar) -t- |iV,0))/V2 can be made orthogonal in a 
two- mode interferometer by a relative phase shift oi ir/N 
§ , and the two- mode interferometer is a special case of 
the multimode interferometer considered here, we know 
that the smallest necessary phase shift (f> is smaller than, 
or equal to, Tr/2N. 

Minimizing (j) above corresponds to an experimental 
situation where we have an arbitrary number of induced 
phase shifts at our disposal (one in each mode). How- 
ever, we easily identify another experimental situation 
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of interest, where we have a given amount of phase shift 
that can be spht and distributed among the modes. If we 
are also allowed to choose the signs of the phase shifts, 
our problem is to minimize ^tot = Sm=i \4'm\- We will 
see that the solution to this problem also gives a tight 
bound for the case where all phase shifts have the same 
sign. 

In accordance with our discussion above, we consider 
an A/-mode interferometer with phase shifts 4)„i — ^m4>i 
where — 1 < Xm < 1, in the different arms. We will also, 
without any loss of generality, number the modes in such 
a way that Ai > A2 > . . . > Am- The two problems for- 
mulated above can then be expressed in a mathematical 
form as follows: We want to find the smallest value of 6 



> > > \ \ 



and 4>tot = Z]m=i \4>m\, respectively, such that 



I*) =0, 



(1) 



for some pure M-mode state \^) with a given average 
energy. In particular, we are, in correspondence with the 
Heisenberg limit, looking for the minimum product of 
our accuracy measure {4> or 0tot) and the average photon 
number {N). 




Fixed 
Mirror 



Fixed 
Mirror 




FIG. 2. A generalized delay-line Michelson interferometer 
for detection of gravitational waves. As the x- and y-arm are 
perpendicular, quadrupole gravitational waves would, under 
certain conditions, induce changes in the mirror separations 
Lx and Ly that are of the same magnitude but with different 
signs [15]. The present setup then becomes equivalent to the 
one sketched in Fig. |l|. Note that we consider a general source 
and detector, whereas the conventional delay-line Michelson 
interferometer used for gravitational wave detection consists 
of a laser source, a beam splitter and a photodetector. 



III. THE ENERGY EIGENSTATES 

In this section, we restrict ourself to the A''-photon 
states. We derive lower bounds on 4> and 0tot, and give 
examples that saturate these bounds. The general case, 
where the states are allowed to be superpositions of dif- 
ferent photon number states, is considered in Sec. IV. 



A. Minimizing (j> with respect to N 



FIG. 1. A sketch of an interferometric setup to probe the 
displacement Aa; of the two-sided mirror at the center. As 
the movable mirror is displaced by Aa; to the right, the phase 
shift in the left mode becomes (jjidt = 27r7?ieftAx/Aocos7, 
where Riaft is the number of times the beam runs between 
the left and center mirror, Ao is the optical wavelength, 7 is 
the angle of incidence, and we have assumed unit refractive 
index. The phase shift in the right mode is given by a similar 
expression, but has opposite sign. In practice, the induced 
phase shift is always limited, e.g., due to the quality of the 
light beam. However, note that we, in principle, can use as 
many modes on each side of the movable mirror as we like. 
By an appropriate choice of the number of reflections and the 
angle of incidence for mode m, we can realize any phase shift 
satisfying \(j>m\ < 4>y where cj> is the greatest possible phase 
shift for a single mode. 



Let \^ n) denote a pure M-mode, iV-photon state and 
introduce the nonnegative parameters = Am — Am- 
The orthogonality condition then becomes 



*Ar) = 0, 



(2) 



where < < 2 and /xi > /i2 > . . . > /xm = 0. Fur- 
thermore, let {|A/^-'^^)}/"i denote the set of all possible 
orthogonal states 



that for each j fulfill 

M 



(3) 
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€{0,1,2,. ..,iV}. (4) 
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(The superscripts are not necessary here, but will be used 
in the next section.) A general pure TV-photon state can 
then be written 



I* 



where J2j=i 

Jn 

E 



N 



„(^)|2 



Jn 



{N) 



(5) 



= 1. Putting this into Eq. (0) gives 
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(AT) 



= 0. 



(6) 



In order for the imaginary part to vanish, the smallest 
definite phase shift we can use is = t: /2N , which im- 



plies that c^^'' = 0, for all j such that X^^^i ^■ 



M 



AN) 



7^0 



or 2N. We must also have = 2 (remember that 
the parameters are labeled in decreasing order). Now, 
let po denote the sum of |c'-^-'p over all j such that 



J2m=i /^™"m'j ~ ^- "^^^ condition for the real part to 
vanish then becomes pq — {1 — po) — 0, i.e., po = 1/2. 
The simplest state that fulfills this has only two nonzero 



coefficients c 



{N) 



An example of such a state is 



|<i>w) = -^(|A^,0,...,0)+e^''|0,0 



(7) 



where i? is an arbitrary real number. This state makes it 
possible to detect a phase shift of tt/2N with certainty. 
As shown above, this is the smallest phase shift that can 
be detected with certainty using TV-photon states. We 
also note that the state (^ can be seen as a represen- 
tation of a two-mode interferometer, since all but two 
modes (arms) are left unexcited. That is, the resolution 
limit (j) ~ 7r/2TV can be reached with a two- mode inter- 
ferometer. 



B. Minimizing (^tot with respect to N 

Now, assume that j^'jv) is an TV-photon state that 
together with the phase shifts 0m — \ra<t> minimize 
0tot = Zlm=i under the condition (Q). Equation (|l|) 
then implies that Eqs. ^ and (^) are satisfied, where 
/i„i0 = <i>m — 4>M- Since the real part of the expression 

(^ must vanish, we have J2m=i Mm^iT] > 0' fo'" some 
j. In order for the imaginary part to vanish, it then fol- 



lows that 4>^m=i t^mn^m] > TT, for some j. Thus, there 
exists an m such that 0m — 0m > t/TV, i.e., 0tot > tt/TV. 
This bound, too, can be saturated by the state (0), if we 
choose 01 = 7r/2TV and 0m = — 7r/2TV. Again, we find 
that the lower bound can be achieved with a two-mode 
interferometer. 



C. The de Broglie wavelength 

Note that the bounds we have obtained for the TV- 
photon states are easily interpreted in terms of their de 
Broglie wavelength AdB = Aq/TV, where Aq is the opti- 
cal wavelength [p| . We know that an ordinary two- mode 
interferometer fed with classical light (or a single pho- 
ton) needs a relative phase shift of tt, corresponding to 
half the optical wavelength, in order to change between 
constructive and destructive interference. If we instead 
feed the interferometer with nonclassical states of light 
with an average energy (TV) , the corresponding de Broglie 
wavelength is limited from below by Ao/(TV). Thus, the 
necessary total phase shift is 0tot — T^/iN). That the 
bound does not depend of the number of modes M, can 
be seen as a consequence of the fact that the de Broglie 
wavelength corresponds to the interfering entities ( "pho- 
ton clusters" ) and not to the number of interfering paths. 
However, one could expect that it would be possible to 
divide this total phase shift equally among all the M 
modes, so that = tt/M{N) would be sufficient. How- 
ever, as we have seen above, this is not true. We found 
that the total phase shift can only be divided between 
two modes, which gives > 7r/2(TV) and makes it possi- 
ble to reach the bound with a two-mode interferometer. 

To simplify comparisons between the limits for differ- 
ent experimental situations, we have gathered our results 
in Table |. The bounds for the single-mode case follow 
from Ref . , and the set of states saturating these was 
derived in Ref. [ |T2| . Since single- mode photon number 
states cannot be made orthogonal by a phase shift, the 
bounds found in this section are the only limits given for 
the TV-photon states in the table. 



IV. THE GENERAL CASE 

Next, we tackle our problem for a general pure state. 
We make some general considerations before we concen- 
trate on the minimization of the respective measures 
and 0tot in two subsections. 

TABLE I. Lower limits on the distinguishable phase shifts. 





Entity 


One mode 


Two or more modes 


Arbitrary states 





Tv/2{N) 


^ 1.38005/(A> 




0tot 


7r/2(A) 


7r/2(A) 


A-photon states 





NA 


7r/2(A) 




0tot 


NA 


7V /{N) 
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Let us start by denoting a general pure state as 

oo 

= ^r^e^^^lM/^), (8) 



7V=0 



N 



where and ipN are real and nonnegative, J2'n=o ' 
1, and the states \^ n) are given by Eq. (||). Our orthog- 
onahty condition ([l]) can then be written 



^ r^crjv = 0, 
where we have defined 



tTjV 



exp Xmul^^j j . 

j=l \ m=l / 

The average photon number of the state (^) is 

oo 



(9) 



(10) 



(11) 



JV=0 



In order to find the normahzed states that minimize the 
energy for a given ip under the condition we define 
the function 



F 



/ . 'N 

N=0 



rlN + a 1 



N=0 



N 



-PYrlaM, (12) 

N=0 



where a and /? are Lagrange's multiphers. Any state 
that minimizes the energy represents an extremum of F . 
Therefore, these states must satisfy 

dF 

— = 2rk{k-a + f3<Jk)^0, Vfc. (13) 
ork 

Since fc, a, and (3 are real, this means that for each k 
such that r/c 0, CTfc must be real. 

Now, let the values of k satisfying rj; ^ be denoted 
kg, and let |e) = I], r;,^ exp(i(^fej|5'fej be an optimal 
state. Then there always exists a largest positive value of 
(1 — CTfc )/kq for those kg that satisfy kg ^ 0. Let iV > 
denote one excited manifold that attains that value, and 
define the state 



|0> = x/l-r2,-<5|0,0,...,0) + 



N 



8W 



where 



6 = 



1 - gfc, 
1-aN 



(14) 



(15) 



It is easily verified that |0) is normalized and satisfies the 
orthogonality condition (^. The average photon number 
is found to be 

N<yrlkg, (16) 



r%N + 



E I 



cfn 



i.e., it is smaller than, or equal to, the average photon 
number of \Q). Hence, all optimal combinations of the 
average photon number and (j), or (/>toti can be attained 
with states that are superpositions of the vacuum and 
one excited manifold N . 



A. Minimizing (j> with respect to (iV) 

In accordance with our findings above, we restrict our 
attention, without loss of generality, to states of the form 
(^), which we rewrite as 



19) =r2|0,0,...,0) + 



l-^ol*Jv) 



(17) 



From the orthogonality condition (^) and Eq. (^l]), we 
then obtain = -rl/{l - r^) = 1 - N/{N). Thus, our 
problem is to minimize crjv, and thereby (A'"), for a given 
phase shift cj). Assuming an optimal state, must be 
real according to Eq. (13). It then follows from Eq. (|l0| ) 
that (Jat > cos(A^0) for phase shifts (p < tt/N. The lower 
bound can be saturated, e.g., by choosing Ai = 1, Am = 
— 1, and the state (|^). The minimum value = — 1 is 
then reached with = w/N. Introducing rj = N/ (N), we 
have 



cos(A(/)) = 1 — 77 



arccos(l — 77) 



(18) 



In order for the phase-shifted and original state to be 
orthogonal, rg < 1/2 must be satisfied. Thus, we have 
1 < ^7 < 2. The value of rj that minimizes ((> for a given 
energy must satisfy 



1 



7]^i - (1 - Tjy 



arccos(l — 77) 



which implies that 



arccos(l — rj) = 



V 



0, 

(19) 
(20) 



Thus, the optimal value of 77 is found to be 

77opt « 1.6891577, (21) 

which gives 

arccos(l — 77opt) 1.38005 



Vopt{N) 



(N) 



(22) 



We notice that the bound can only be saturated for 
(A) ~ N/T]opt, where A^ is a positive integer. This can 
be achieved by choosing, e.g., Ai = 1, Am = — 1, and the 
state 



' '7opt - 1 



Vopt 

1 



|0,0,...,0) 



(|A^,0,...,0)+e^''|0,0,...,A^)) 



\/2?fcpt 

« 0.638740 |0, 0, . . . , 0) + 0.769423 |$7v), (23) 

where 7? is an arbitrary real number and |$jv) is defined 
in Eq. (^. Again, we see that two modes are sufficient 
to obtain the best possible accuracy. 
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B. Minimizing ^tot with respect to (N) 

The minimum value of 0tot for a general state of a 
given average energy follows directly from the work by 
Margolus and Levitin ||ll| . They derived the time of free 
evolution necessary for a state to evolve into an orthog- 
onal state. For the optical field, the necessary time is 
T = 1/4/ (iV), where / is the optical frequency. Since 
this corresponds to a phase shift — tt/2{N), we must 
have (/)tot > n/2{N). The bound can be saturated, e.g., 
by using the state 



|r!^) = -i=(|o,o,. 



,0) + |iV,0,...,0)) 



(24) 



and choosing (pi — tt / N as the only phase shift. We see 
that this bound, too, can be achieved with a two-mode 
interferometer. In fact, it can even be achieved with a 
single mode. For further comparisons, we refer to Table |. 



V. ENERGY MEASUREMENTS AND 
WELL-DEFINED RELATIVE PHASE 

In the previous sections, we have derived the small- 
est distinguishable phase shifts under two different re- 
straints. We know from quantum theory that each time 
we make a measurement corresponding to a Hermitian 
operator with the original and phase-shifted state as 
eigenstates, the outcome will tell us whether the phase 
shift was applied or not. However, it is, in general, hard 
to realize such a measurement. In fact, most often the 
energy is measured. For example, in quantum optics it is 
usually the photon number, or the intensity of the opti- 
cal field, that is finally probed. In Fig. |^, we have drawn 
a general measurement setup of this kind. Such mea- 
surements would, in principle, allow us to achieve the 
resolution derived for the energy eigenstates in Sec. [O. 
In contrast, the vacuum component of the state (E3) is 
invariant under phase shifts, and makes it impossible to 
tell if the state has been phase-shifted or not when no 
energy is found in the measurement. 

Ideally, the interference, which is represented by 
Kietector in Fig. ^ is lossless and unitary. In a mea- 
surement that gives the total energy, the requirement to 
be able to distinguish a phase shift (p with certainty can 
therefore be written 



(*|ljve 



l*>-o. 



where we have used the definition 



1 



N 



Jn 



(25) 



(26) 



That is, for any possible outcome of the total energy, the 
original and phase-shifted states must be orthogonal. 
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FIG. 3. A schematic diagram of an interferometric mea- 
surement that also gives the total energy of the measured 
state. The interference of the M modes is characterized by 
the lossless and unitary transformation t/dctcctor- 

In an earlier paper [ p^ , we claimed that if |^) is a two- 
mode state, /7ps(0) = exp(i0[fii — fi2]/2) is the two-mode 
unitary relative-phase shifting operator, and if 



(eic/ps(0)io -0 



(27) 



can be fulfilled for some value of </>, then "Eq. ( P7|) consti- 
tutes the mathematical criterion for a state with a well- 
defined relative phase." It is clear that Eq. ( p7| ) is a spe- 
cial case of Eq. (^ , and therefore constitutes the neces- 
sary and sufficient condition for a relative phase shift to 
be distinguishable. However, as the following example 
shows, this property should not be referred to as "well- 
defined relative phase." Suppose that 



\0 



1 

71 



(|0,0) + |0,1)) 



1 



|0)®^(|0) 



|1)). (28) 



It is trivial to show that in this case Eq. (|2^) is fulfilled for 
(j) = 27r. However, a product state where one of the con- 
stituent factors is a single mode vacuum state is hardly a 
sensible candidate for a state with a well-defined relative 
phase. This is confirmed by the relative-phase distri- 
bution functions, which are flat |17|. The conclusion is 
that the criterion for an operationally well-defined rela- 
tive phase must be reformulated. 

The constant relative-phase distribution functions of 
the state ( p8| ) are a consequence of the fact that the 
relative-phase and the total energy are compatible ob- 
servables [Q. Since Ups{4>) also commutes with the total 
photon number operator X]w=o ^^n, any measurement 
of two-mode relative phase (or visibility) will reveal the 
total photon number. Therefore, in order for the phase 
shift (p to be perfectly distinguishable, the phase-shifted 
state C/ps(0)IO niust, as in Eq. (p5|), become orthogonal 
to the original state |^) in every photon- number mani- 
fold. Consequently, the definition of well-defined relative 
phase should be refined as follows. If, for some phase 
shift (j) and for all N, it is possible to fulfill 



(29) 



then the state |^) is said to have a well-defined relative 
phase. That any state fulfilling Eq. (E9|) also fulfills Eq. 
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( P7| ) can easily be verified. However, the converse is not 
true, as shown by the state (^8|), which does not fulfill 

Since all the components in different excitation mani- 
folds must b e ma de orthogonal simultaneously, the limits 
found in Sec. |lll| imply that in this case (j> > 7r/2iVmax and 
0tot > Ti'/^max, whcrc iVmax dcuotes the highest excited 
manifold. 
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VI. DISCUSSION 

We have found tight bounds set by the average energy 
of the used states on our accuracy measures (j) and (/)tot ■ 
The bounds were derived for optical interferometers and 
pure states, but are obviously valid for other interfering 
bosons and mixed states, too. Our results also give tight 
bounds in the case where all phase shifts have the same 
sign. This would correspond to an experimental situation 
where there is a limited amount of phase-shift inducing 
material that can be distributed among the modes. Since 
the state can saturate the bounds for (jjtot with non- 
negative or nonpositive phase shifts (jjm alone, the magni- 
tude of the sum of these phase shifts has the same bounds 
for TV-photon states and general states as 0tot- 

We would also like to point out that our results were 
derived from fundamental considerations. Apart from 
the brief discussion in Sec. 0, we did not consider the 
problem of how to realize the measurements necessary to 
distinguish the optimal states. Although we have found 
that there is no fundamental advantage in using mul- 
timode interferometers in terms of their accuracy, they 
may offer simpler realizations compared to two-mode in- 
terferometers. For example, it may be easier to generate 
multimode states that give the same resolution as the 
two-mode states presented here. 

Let us also comment on the earlier investigation of 
multimode interferometers by D'Ariano and Paris [ p^ , 
which arrived at the conclusion that the multimode in- 
terferometers actually improve the accuracy. They con- 
sidered interferometers whose phase shifts in the differ- 
ent modes were 4>m — nriLp. The "variance" of the es- 
timate of ip was found to be inversely proportional to 
the number of modes A/, indicating that the accuracy 
in the measurement of was improved correspondingly. 
However, we note that the largest relative phase shift 
4'M — 4'i = — l)<p increases with the number of modes 
M. Since the two-mode interferometer employs a rela- 
tive phase shift of (p instead of (M — l)(p, its sensitivity 
is not maximized. The relative phase shift {M — l)ip 
would obviously give the same improvement in the accu- 
racy of if for a two- mode interferometer as D'Ariano and 
Paris found for the Af-mode interferometer. Therefore, 
we conclude that for a given average energy the ultimate 
phase-shift resolution of an interferometer is independent 
of the number of interfering modes. 
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